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We demonstrate that the time operator that measures the time of arrival of a quantum particle into chosen 
state can be defined as a self-adjoint quantum-mechanical operator using periodic boundary conditions on ap- 
plied to wavefuncions in energy representation. The time becomes quantized into discreet eigenvalues and the 
eigenstates of the time operator, the stroboscopic wavepackets introduced recently [Phys. Rev. Lett. 101, 
046402 (2008).] form orthogonal system of states. The formalism provides simple physical interpretation of the 
time-measurement process and direct construction of normalized, positive definite probability distribution for 
the quantized values of the arrival time. The average value of the time is equal to the phase time but in general 
depends on the choise of zero time eigenstate, whereas the uncertainity of the average is related to the traversal 
time and is independent of this choise. The general fromalism is applied to a particle tunneling through resonant 
tunneling barrier in ID. 

PACS numbers: 03.65.Xp, 73.63.-b, 72.10.Bg 



I. INTRODUCTION 

The concept of time operator in quantum mechanics is a 
difficult and confusing one|[l}[3l. Heisenberg formulated the 
time-energy uncertainity principle already in the early days of 
quantum theory, indicating through analogy with the position- 
momentum uncertainty principle, that some time operator 
should exist. However, shortly after this, Pauli argued that no 
such self-adjoint operator can exist 0. Further development 
in scattering theory pursued the search for an estimates of 
time-scales associated with quantum processes, establishing 
the phase time, the time delay seen in the motion of the max- 
imum of a wavepacket as the relevant quantity fS', "Fl. This, 
however, turned out to be unsatisfactory due to inherent ambi- 
guity in the preparation of the wavepackets or identification of 
its maxima or other features. Several imaginative approaches, 
like the so called Larmor-clock time ^ IT HTOl or the traver- 
sal time 111 [TT] |l2| were suggested to identify the relevant 
time-scales. However, no final formulation of the problem 
has been established nor a consensus has been reached if such 
a formalism should exist. Nonetheless, the time-scale related 
to tunneling is an extremely useful concept for relevance of 
many-body effects in electronic transport through nanostruc- 
tures. This has been analyzed in the pioneering work by Jon- 
son 1 13 , 14| where the time scale of tunneling was determined 
by the time scale of formation of the image charge, causing 
alteration of the effective tunneling baiTier ifTSl [T6l . More 
generally, our ability to characterize the time-scales of tran- 
sit or tunneling of an electron through a nanocontacts would 
be extremely helpful in understanding the importance of inter- 
actions in ab initio description of quantum transport ifTTl - ISSl . 

Independently of these physically motivated treatments, a 
important step forward in understanding the time operator, 
not as a self-adjoint operator but rather as a positively valued 
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operator measure, has been done by Holevo fZA]. Indepen- 
dently, Kijowski 1.25,1 heuristically constructed a distribution 
of a time-of-arrival for a quantum particle. His work was later 
developed into the formulation of the construction of probabil- 
ity distribution of the arrival time|26 - 29 1. Most recently, even 
the problem of non self-adjointness of the time operator for a 
free quantum particle has been addressed by Galapon f30l l3TI 
by introducing the confined time of arrival operator (CTAO) 
within finite space using specific boundary conditions in the 
real space. 

In the present paper, we propose an alternative formalism, 
the use of periodic boundary conditions in the energy repre- 
sentation that leads to a family of self-adjoint time operators. 
In contrast to the CTAO, arbitrary scattering potentials can be 
considered from the start and the related issue of normaliza- 
tion of the probability distribution for times 1261 l27l is also 
resolved. The boundary conditions in the energy representa- 
tion lead to formal quantization of the time, similarly to the 
situation with the CTAO. The quantization of time is a use- 
ful mathematical tool for drawing a simple physical picture of 
the time-dynamics of the quantum paiticle within the orthog- 
onal time-eigenstate basis. This is used for physical interpre- 
tation of the zero-time eigenstate and its relation to the con- 
ventional arrival-time operator [3| and the 'time-of -presence' 
operator ll24l . However, within the energy representation, 
the formalism is similar to many previous uses of the time- 
operators in the form of diferentiation by energy |6, 24, 32]- 
l35l . Our formalism is demonstrated on a simple example of 
scattering of a particle on a resonant potential in ID. 



II. DEFINITION AND GENERAL PROPARTIES OF THE 
TIME OPERATOR 

In our work we will consider a quantum particle moving 
along X axis, characterized by its Hamiltonian H assumed to 
have continuous spectrum, occupying a state | ) . For this par- 
ticle we introduce a family of time operators, f,;, where each 
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of them gives the time it took for the particle to arrive into the 
state \ assuming its dynamics has been governed by H. The 
concept of the time operator also demands a definition of the 
zero of the time for certain state, which will be discussed af- 
ter the operator is introduced. Different zero-time states are in 
one-to-one correspondence to different members of the family 
of time operators. 

For the state we will assume that it can be expressed as 
a linear combination of the Hamiltonian eigenstates, \e,a), 
with energies from a finite interval, the energy band e E 
(eo,eo + Ae), 



I0) = L 



£o+Ae 



(1) 



where a is the quantum number for degenerate states at the 
energy £. This degeneracy arises from some other operator A 
(or possibly several operators) that commutes with the Hamil- 
tonian. The eigenstates are chosen so as to be common eigen- 
states of both A and H. The complex amplitude (e, a) rep- 
resents the state |0) in the energy-band representation. The 
eigenstates are normalized to the delta function of energy, and 
the amplitude is normalized to one. This representation is not 
unique, another one can be obtained choosing a different op- 
erator B that also commutes with H (but not with A). The 
representations are then related through an energy-dependent 
unitary transformation 



|e,a)=^f/„,p(e)|e,j3). 



(2) 



We will assume that the amplitude for the state |0) has a 
support within the considered energy band so that at the ends 
of the interval we have g^ (go, a) = g,;, (eq + Ae, a) =0. The 
states g^{e,a) form a subspace of the Hilbert space of 
all square-integrable states g{e,a) that are periodic within the 
energy band. Extending the width of the energy band, or con- 
sidering a union of all energy bands covering the continuous 
spectrum of the Hamiltonian, and adding its possible bound 
states forms a complete set of states ll36l . However, for the 
purpose of introducing the time operators for state \(j)) it is 
sufficient to consider single energy band. 

We will demonstrate that in the energy-band representa- 
tion, within the self-adjoint time operator can be defined 



de 



(3) 



Apart from the so-far unspecified Hermitian, energy- 
dependent matrix TJa.a'l^)' ™d the fact that we define it only 
within the energy band, this operator has been known for long 
time as the operator for time in the energy-representation. It 
is well known that it is not self-adjoint if the whole spec- 
trum is considered |24| and not unique by the freedom of 
choice in the energy representation ll28l . The former is re- 
moved by the finite energy interval and the periodic bound- 
ary conditions employed. On the other hand, the freedom of 
choice of the energy representation in Eq. [T] is related to the 
choice of the Hermitian matrix in the definition in 



Eq. [3] Using the unitary transformation, introduced in Eq. |2] 

l^a.pi^) — 6xp{//'^(ie'Va^/3(e')} we find a transformed time 
operator 



iSa,a' ^ + J7a,a'(e) " Va.a' (e)- 

de 



(4) 



Hence, starting from a particular energy-representation and 
a particular choice of J]^, we can find a representation 
where the time operator is represented by the energy derivative 
only. This latter representation, if we had some rationale for 
choosing it independently of the time operator, could serve 
as the basis for definition of the time operator without any 
ambiguities. 

The first step along this line is to demand that the time op- 
erator should commute with chosen operator(s) A. Examples 
of these could be the linear or angular momentum, spin, etc. 
One of these is also the projector to the right- and left- going 
scattering states leading to a — R ot a ^ L in ID scattering 
that will be used in the next section. This reduces the matrix 
rja,a' (s) into diagonal form and specifies the time operator for 
a processes which conserve the particular quantum number a. 
For the simplicity of notation, we will not indicate this fact 
with any additional index for the time operator f,j . 

Further specification of fjj is related to the choice of zero 
time state, as discussed below, but in general no unique defini- 
tion of the time operator will be given. Instead, we will accept 
that we deal with a family of operators of the form given by 
Eq.|3]with J7a.a'(e) = 5a.a'^a(£) and that for a specific cal- 
culations we need to choose one particular form. 

The eigenfunctions of the time operator are 



where 



2n 

T„, = — m, m = ± 1 , 
Ae 



(5) 



(6) 



are discreet eigenvalues of the time operator. Similarly to our 
finding, discreet eigenvalues were found in the construction of 
the confined time operator |30 | for a free quantum particle. In 
both cases, the quantization is nothing fundamental and arises 
only as a result of the choice of the energy band: the size 
of the time-quanta could be changed by simply changing the 
width of the energy band while the final average values of the 
time operator will remain independent of this choice, which 
becomes obvious when using the energy representation. Still, 
choosing a wider energy band, i.e. decreasing the quantum 
of time, is desirable if one is interested in finer details of the 
probability distribution in time variable. 

Rewriting the time eigenstates from the energy-band repre- 
sentation into the abstract form we have 



T,„,a = 



1 



Ae Jeq 



t/ee-'"'"V'"''^'''«('=')|e,a) 



(7) 



This set of states is subspace of the stroboscopic wavepacket 
basis, recently introduced for the description of open non- 
equilibrium electronic systems 1361 l37l . Here we see that it 



3 



naturally arises as the set of eigenstates of the time operator 
defined on the chosen interval of energies. 

The eigenstates with m — Q, |to = 0,o:), have eigenvalue 
of the time operator zero, i.e. this is the choice of the zero 
of time. Due to the unspecified phase J7a(e) there is a certain 
freedom in the choice of this zero time eigenstate (and hence a 
particular time operator). If a particle is in one of these states, 
it took zero time to arrive into it. On the other hand, for a 
state |T„,,a) it took it precisely the time T,„ — 27Zm/{Ae) for 
the particle to arrive there from |to,o:), since from Eq. |7]we 
find 



-iHr,,, 



10, a) = \T,„,a) 



(8) 



Finally, for a particle in an arbitrary state within the energy 
band, |0) = Lm,a Cm,a Km, a), and |cm,aP= |(Tm,a|0) p will 
be the probability for the particle to arrive into it in time = 
2nm/{Ae). 

From the above it follows that the expectation value of the 
time xfi in the state |0) is given by weighting the different time 
eigenvalues with the probabilities that the relevant eigenstate 
is present in the state 



4 = L km,apT,„ = (0|f^ 10), 



(9) 



This is equivalent to using the form in Eq. [3]within the energy- 
band representation. 



de 



i(e,a). (10) 



which motivates the formal definition of the time operator by 
the Eq. [3] Clearly, the expectation value of the time opera- 
tor depends on the choice of zero-time eigenstate, i.e. on the 
choice of the phases J7a(e). In contrast, for the uncertainty of 
this average. 



we find 



(11) 



(12) 



which is manifestly independent of the choice of the phases 
77a(e) and hence characteristic of the whole family of time 
operators. 

The whole family of time operators, Eq. [3] fulfills the 
canonical commutation relation with the Hamiltonian H = 
„/, [f,H] = /, if the latter is understood to act only on 
states 10) with finite support within the energy band. (A mi- 
nor technical issue that can be dealt with arises if the whole 
i^ti£ is considered, since there the Hamiltonian is not continu- 
ous at the boundaries of the energy band.) This commutation 
relation then leads automatically to the uncertainty relation 
for the mean square fluctuations in the energy and the time, 
AitsH > 1 /2. 

The argument due to Pauli |3 4| regarding the non- 
existence of the self-adjoint time-operator does not apply 



since the boundary conditions cause the energy-shift operator 
to move the states periodically within the band. Namely, using 
the orthogonality of the time-operators' eigenstates (Eq.[5|) we 
can expand the Hamiltonian's eigenstates \£,a), Eq. iTTinto 
the former and find the identity 

He-'^'' |e, a) = [(£ - £')mod(A£)] e"'^'* |e, a) , (13) 

for £ e (£o , £o + A£ ) . On the other hand, if the periodic bound- 
ary conditions within the bands were not used, the above iden- 
tity would not contain the modulo operation with the differ- 
ence (£ — £') and the result would be that the state e^'^ ^ |e, a) 
is an eigenstate of Hamiltonian with the eigenvalue £ — e'. 
Following Pauli, and in view of arbitrariness of £', this would 
be in contradiction with the existence of the lower bound on 
the eigenenergies. However, we have shown above that the 
use of the periodic boundary conditions removes this prob- 
lem. 

The here used time operator is, by its character, close to 
the 'time-of-presence' mentioned in the review by Muga and 
Leavens |3|. However, many authors IZSHZTl l29l prefer the 
concept of the arrival-time operator that gives the average 
value of time for a quantum particle to arrive at a spatial po- 
sition XQ if initially (at time f = 0) it was in a chosen state 
I xjf) . One can easily see that such an operator is given by — f,; , 
with a specific choice of the phases J7a(£). The latter is such 
that the zero-time eigenstates resembles the position eigen- 
state 8{x — xq) as much as possible. For example, for a free 
quantum particle the energy eigenstates are 



U £,a = 



1 



;t= V2£ 



(14) 



and using the projection /t/£ |£, a) (£, a|xo), one finds 
^a(e) = {d/de)y/2exo. The interpretation of this time op- 
erator is as follows: we expand the state |i//) into the arrival- 
time operator eigenstates, jy/) = {%n;C^\w) K;n,0!)- Then 
I (T„,,a|v/) 1^ is the probability that the particle in jy/) will ar- 
rive into the |T() ~ 0, a) in time — T„,. Identifying the zero-time 
eigenstate with measurement device at x = jcq gives the sought 
Kijowski probability distribution Il26ll27l . However, we need 
to stress that the arrival state, i.e. the zero-time eigenstate, 
can be quite different from the position eigenstate 5{x — X()) 
so that it should not be interpreted literally as the probability 
of the time of arrival into xq exactly. 



III. ARRIVAL TIME IN TUNNELING PROBLEMS 

We will now demonstrate the use of the time operator for 
calculation of the tunneling time scales involved in the dy- 
namics of quantum particle in ID. For the state into which we 
expect the particle to arrive we initially take a state |0), lo- 
cated on the right of the tunneling barrier and characterized 
by a momentum directed away from the barrier. 



(x|0)=y^ -^A^{£)e''<^^-^^\ Xfi»0, (15) 
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FIG. 1: (color online) The form of the potential energy used for 
demonstration of the tunneling time scales in ID. The potential has 
two delta-shaped barriers of strength A at its both ends and a constant 
value u in between. Varying the latter gives access to various trans- 
port regimes - from resonant tunneling to opaque tunneling. The 
energy-band of the stroboscopic wavepacket representation is indi- 
cated with arrows. 



where xr determines the average position of a particle in the 
state and k — y/2e. The real amplitude A|ji(e) is a continu- 
ous, differentiable function with its support within the energy 
band. The time it takes for the particle to arrive into the state 
10) from the left of the banier will contain contribution of the 
time it took for the particle to tunnel. 

The dynamics of the particle is governed by the Hamilto- 
nian H which asymptotically, for x — > ±°o, is that of a free 
particle. Close to the origin there is non-zero potential en- 
ergy V{x). The Hamiltonian posses a continuous spectrum 
of doubly-degenerate energy-normalized right- and left- going 
eigenstates {x\\ifg K) and {x\\i/gx) of which we will explicitely 
need only the right-going ones. 



{x\W£,r) = WeAx) 




(16) 



where k — v2e and r(e) and t{e) are the reflexion and trans- 
mission amplitudes respectively. 

As a illustrative example, we will consider the Hamiltonian 
H = {l/2)d^/dx^+V{x), where V{x) X{5{x) + 8{x-a)) + 
(m/2)(1(x) — l{x — a)) where l{x) is the unit-step function, 
X = 1, a = 10.0, and v = 0.1 — 0.65 is a variable potential 
within the delta-functions (see Fig [T]). As can be inferred 
from the amplitude and the phase of transmission amplitude 
shown in Fig. |2] these values offer variety of different trans- 
port regimes. It might be also interesting to mention that this 
Hamiltonian coiTesponds to a simple model of a perturbed 
monoatomic sodium chain ll37l . For the arrival state (j){x) 
in Eq. 15 we take A(e) — Na [cos^ ((e — ei);r/e2)] , where 



ei = eo + Ae/2 and £2 = Ae, with eo = 0.2 and Ae = 0.4, 
which produces a convenient localized state covering the in- 
teresting features in the transmission and its phase (Fig. |2|l. 
The constant Na is fixed by the normalization of the state (p (x). 
The localization of the arrival state will be set at xr = 100 
which for the chosen parameters will guarantee that its ampli- 
tude in the region of the nonzero potential is negligible. 

We proceed by the selection of the time-operator Firstly, 
the Hermitian matrix T7a,a'(e) — J?a(e)5a.a' shall be diagonal 
in the basis of the scattering states so that a=Roi L for right- 
or left- going scattering states. Secondly, the phases of this 
matrix will form the zero-time eigenstate from the incoming 
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FIG. 2: (color online) The amplitude and the phase of the transmis- 
sion amplitude f(e) for a scattering state in the potential given in 
Fig. [TJfor u = 0.3 (upper panel) and u = 0.65 (lower panel). Since 
the energy band consists of energies e e (0.2,0.6), the smaller u cor- 
responds to resonant transport where as the larger ;/ gives dominantly 
opaque tunneling regime. 



scattering states at some initial time Tq, 
lR/Li£) = -To- 



(17) 



If To were << 0, the zero-time eigenstates would look just 
like wavepackets formed from incoming plane-waves, local- 
ized far to the left (right) of the barrier for a — R (L) respec- 
tively. However, this only shifts its origin for the time, so that 
in our calculations we will simply use Tq = 0. 

To evaluate the average time of arrival into |0), we need 
to express this state in the energy representation of the above 
scattering states, 

\(P)=Y, fde\xi/e.a){VfeM<t>), a ^ R,L (18) 



For the chosen state \<j)), Eq. 15 both 
will be nonzero 



{We.r\<P) and (v/e.il^) 



The presence of the left-going states goes 
against our intention to characterize the tunneling time scale. 
The physically relevant state in which the time average should 
be calculated should consists of the right-going states only. To 
construct the final state correctly for a particle moving from 
the left to the right of the barrier, we need to construct a state 
1 0), obtained from the state |0) by a von Neumann projection 
on the right-going states only, / de \ \I/£,r) {Y£,r\' 



1 



de\\(f£A{\lf£,R\<j)), 



where is the normalization constant. 



N= / de\{We.R\<P)\ 



(19) 



(20) 



The state for = 100 is negligibly small in the region 
where the states yfR{x) differ from their asymptotic form for 
.X >> in Eq. 16 so that 



(V/£,r|0) =f*(e)A^(e)6 



(21) 
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FIG. 3: (color online) Probability densities of projected states 0(x) 
for x« = 100 into which the particle arrives, for two different poten- 
tial barriers u. In contrast with the unprojected state (p{x), the arrival 
states are distorted and contain significant weight on both sides off 
the potential barrier (the amplitude for the state for ;( = 0.65 on the 



right is 10^ 
xe (0,10). 



times magnified so that it is visible) localized within 



andA^ = / de\A{e)t{e)\^ = (T) has the meaning of the average 
transmission probabiHty of the particle from the left to the 
right. The unprojected (j){x) and the projected states 0(x) for 
two extremal values of the potential parameters u are shown 
in Fig. |3] It should be noted, that unlike the projected 
states are distorted and contain a non-zero amplitude of the 
reflected state on the left of the barrier This is necessary to 
have the total probability of arriving into the projected state at 
any time being one. 

According to the general treatment, the probability that the 
particle arrived into state 1 ) at time T„, is 



P,„ = \{t,„,R\^)\^ 



1 

Jf) 



de' 



(22) 



which is positive definite and normalized to one. The example 
of several such probabilities obtained for our particular exam- 
ple is shown in Fig.|4] For the energy band above the potential 
between the barriers (m = 0. 1 , 0.3), the particle will bounce be- 
tween the barriers, resulting in the probabiUty density of the 
arrival time into state to the right of the potential with several 
local maxima, separated by a time it takes for the particle to 
traverse the distance between the barrier twice. On the other 
hand, for a particle within a energy band below the poten- 
tial barrier, no particular structure is visible. Clearly, different 
potentials lead to different distributions for arrival times and 
only the full probability distribution gives a complete picture 
of the time scales involved in the dynamics. Still for many 
cases two measures are most useful - the average value and 
its uncertainty. 

The average time of the particle to arrive into the state 1 0) is 
most easily calculated within the energy-band representation. 



1 



de\A^{e)t{E)\^ 



dO xr 

■1- + — 
de v„ 



7b 



, (23) 
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FIG. 4: (color online) The probabilities of arrival into state ^(x) 
at time Tm for different transport regimes. In resonant transport 
{u = 0.3), due to good localization of the state we see several max- 
ima giving multiple bounces between the delta-function barriers; in- 
creasing u into tunneling regime we first see very broad distribu- 
tion (u = 0.53) with its width characterized by the Buttiker-Landauer 
time, evolving into relatively narrow saturated distribution arriving 
earlier than the free particle (the Hartman effect). 



where 0(e) is the phase of the transmission amplitude t{e) ~ 
|f (e)|e'^. The additive term proportional to the average posi- 
tion Xr is the classical expression xr / Vg, where Vg ~ de/dk is 
the group velocity of the particle outside of the barrier. The 
first term is the generalization of the phase time since for a 
state within a narrow in energy band (Ae 0), the average 
time equals the well known expression dO{e)/de. The choice 
of the phases ?]«(£) = — Tq leads to a simple shift in the time, 
independent of the transmission or the Hamiltonian's poten- 
tial. Keeping this form one can compare average times for 
different scattering potentials V{x), localized close to the ori- 
gin. In general, however, the energy-dependent phase leads 
to nontrivial change in the average time. Simple, potential- 
independent shift is found only if this energy dependence is 
negligible, i.e. (d/d£)riR{e)Ae < Tjff(e). 

According to Eq. [12] the uncertainty of the average time is 
given by 



iATf = J^{\A,ie)tie)\'f. 



(24) 



There are two limiting cases: (1) {d/de)A(^ 3> {d/de)\t\ the 
uncertainty is dominated by the energy width of the state | 
A^(e), and it does not carry information about the time scale 
in the scattering, and (2) {d/de)A^ <C {d/de)\t \ when the un- 
certainty is dominated by the time scale known previously as 
the traversal time [II \ and identified [TJ as one of the Lar- 
mor clock times, = |f |^'((i/iie)|f |. Most importantly, the 
traversal time dominates the uncertainty whenever the the en- 
ergy of the state moves under the barrier and the transmission 
is principally given by \t\ ~ exp{—d^/2{u — e)}, where d is 
the barrier width and u its heigh. In Russian literatures it has 
been known as the Keldysh time| 12 , 38 1 and is approximately 
given as Az ^ d/K where K — \J'l{u — e) is the magnitude 
of the imaginary momentum under the barrier Our result are 
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also in agreement with the observations of Yucel and Andrei 
that this is the time scale that should appear in the energy-time 
uncertainty principle ||38l . The identification of the traversal 
time with the variance of probability density of the dwell time 
has been also found by Olkhovsky [32J . 

The appearance of the huge uncertainty in time, At, due to 
presence of the average energy of the state just below the bar- 
rier is demonstrated in Fig. |4]for u = 0.53, in the very long 
tail towards large times of arrival. This will have significant 
effect on the tunneling process itself if the particle could inter- 
act with some additional degree of freedom within the barrier 
This explains the appearance of this time scale in interacting 
tunneling models yTl[T3lfT4ll38]| . On the other hand, increas- 
ing the potential barrier further, the uncertainty is dominated 
by the energy width of the state ^{x), and hence independent 
of u (e.g. u ^ 0.55 and u ~ 0.65). 

Similarly to the traversal time, one can find a close cor- 
respondence between the average of and the Buttiker- 
Landauer time 1 1 1 Tj-, but this statement, as any expression 
involving the phase time, is valid only for a particular choice 
of the phases in the zero-time eigenstate. This, however, is 
less important in tunneling regime, where At is the dominant 
contribution to the average value of . 

In the conclusions, we have introduced a family of self- 
adjoint time operators into the framework of standard quan- 
tum mechanics using periodic boundary conditions used on 
the amplitudes within the energy-band representation. This 



representation leads to quantization of time which is a use- 
ful tool to regularize the time-eigenstates and to draw physi- 
cal interpretation of the use of this operator. We have shown 
that each member of the family of time operators fulfills the 
canonical commutation relation with the particle's Hamilto- 
nian and hence the energy-time uncertainty principle, identi- 
fied its eigenstates are stroboscopic wavepackets, and showed 
how the present treatment avoids Pauli's argument. In the con- 
text of the tunneling of a quantum particle, we have shown 
how a positively defined distribution function of times of ar- 
rival into an arbitrary state is obtained. Using the latter we 
have obtained the average value of the time operator which 
corresponds to the phase time which, however, is specific to 
a particular choice of the zero-time eigenstate and hence the 
particular choice of the time operator. In contrast, its uncer- 
tainty is independent of this choice and in the limit of narrow 
energy spread of the state it is equal to the traversal time scale 
of Buttiker and Landauer, or the Keldysh time. Our formalism 
confirms the role of energy derivative in energy representa- 
tion as a legitimate time-operator in non-relativistic quantum 
mechanics and opens consistent ways for studding temporal 
behavior in many quantum-mechanical problems of interest. 
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